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NOMENCLATURE 


b temperature coefficient of viscosity 
Cp specific heat of the fluid at constant pressure 
D barrel inner diameter 
G Griffith number = p, V^^/[Kf{Ti, — T,)] 

H height of the screw channel 

Hj, depth (below the screw surface) at which the coordinate system is fixed 
K f thermal conductivity of the fluid 
Kg thermal conductivity of the screw material 
L axial screw length(metering section) 

N screw speed, r.p.m. 

n power law index 

p pressure 

p reference pressure [ = p (Vi^/H)] 

Pav average pressure in the x-y plane at each z-location 
Pe Peclet number (= U^H/o:) 

Q total volumetric flow rate 

dimensionless volumetric flow rate or throughput parameter [= Q/HWHz] 
T temperature 

T, temperature at the inlet of the heating zone 
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Greek Symbols 


reference temperature 
velocity component in x direction 
tangential velocity of the barrel (= ttDN / 60) 
component of Vt along x [= 14sin(^] 
component of Vb along z [= VbCos(f\ 
velocity component in y direction, 
width of screw channel, 
velocity component in z direction 
coordinate axis normal to screw flight, 
coordinate axis normal to screw root, 
coordinate axis along the screw channel. 


aj 

7 

7* 

io 

0 

0bulk 


fio 

P 

<f> 


thermal diflfusivity of the fluid = [ KjIpCp] 
non-dimensional b [ = b(rb-T’,)] 
strain rate 

[= + (g)'} + (g + gf + (gf + 

dimensionless strain rate [= 'yH/Viz] 
reference strain rate 

dimensionless temperature [ = (T - Tt)l{Tb - T,)] 
dimensionless local bulk or average fluid temperature 

reference viscosity 

average viscosity as defined in Eqs. (2.8), (2.10) and (2.12) 
non-dimensional viscosity = pjfi 
density of fluid 
screw helix angle 
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Subscript 
b barrel 

dev developed 

f fluid 

i inlet 

o reference quantity 

ref reference quantity 

s screw 

Superscript 

* dimensionless quantity 
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Chapter 1 


INTRODUCTION 

1.1 Introduction 

The process of extrusion consists of converting a suitable raw material into a product 
of specific cross-section by foicing the material through an orifice or die under 
controlled conditions. Screw extrusion, one such process, is used in many industries 
such as those related to plastics, polymers, pharmaceuticals and food products. In 
order that this conception has practical value, some requirements must be satisfied 
both as legaids the equipment and tlie law material. The equijirncnt must be capable 
of pioviding suitable |)ies.sur(', < ontinuously and uniformly, on the material so that 
it attains icviuin'd pi o|)erti('s. 'I'lins it is necc'ssary to understand the underlying 
mechanism and determine heat transfer rates, residence time distribution and the 
thermal and fluid flow fields etc. within the extruder. 

The screw extrusion machine (Fig. 1.1) consists of a screw of special form 
rotating in a closely fitted heated barrel or cylinder in which feed port is placed 
radially or tangentially at one end and an orifice or die axially at the other. It 
usually has three different geometric sections: a feed section with. relatively deep 
channel depth or flights with greater than normal pitch, a tapered compression or 
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transition section - where temperature and pressure are raised, and a relatively shal- 
low metering section of constant depth - where temperature and pressure increases 
further and high pressure enhances mixing. 

1.2 Fundamental Differences between Polymer and 
Food Processing 

Although there are quite a few similarities between food and polymer extrusion, 
there are some points in which they distinctly differ. They are as follows: 

1. A polymer (or a thermoplastic) can be melted and then on solidification the 
original material can be restored Thus the process is reversible. On the other 
hand, there are irreversible cooking reactions in the food processing. 

2. 'I'lic cooking of food inside an extruder is a complicated process. Both chemical 
and physical changes may occur at the same time under the influence of moist 
heat, pressure and shear. Two most important reactions in food extrusion are 
protein denaturation and starch gelatinisation. 

3. The dough viscosity greatly increases as the reaction proceeds. In contrast, the 
apparent viscosity of thermoplastic will usually decrease as the temperature 
increases. 

RomscMi and Claik (1978) found that for a 25% soy flour suspension, the apparent 
viscosity (.lid decrease initially upon heating until around 1()0"F(71.11°C) and then 
the viscosity started to inciease because of denaturation reaction (here the starch 
gelatinisation reaction is unimportant) demonstrating the cooking phenomenon. 
Thus, during cooking the dough will first encounter “melting” with, a decrease in 
viscosity until a critical temperature is reached at certain point down the channel 
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when its apparent viscosity will show a rise. In short, the viscosity of a food system 
is not only a function of shear rate and temperature, but also depends on the 
composition and the time-temperature history of the process. 

I'hus it is in the terms “melting” and “cooking” that the analogy between a food 
and plasticating extiudcr stops. 


1.3 Literature Review 

Although extrusion as a manufacturing process for structural material can be con- 
sidered to have originated at the end of 19th century, it is only within the last 50 
ycais or so its full potential has begun to be realised. First reported industrial 
machine for extrusion of lead pipe was constructed by John Bramah in 1795. From 
then until the invention of specially designed thermoplastic extruder in 1953 by Paul 
'I'roester in Geimany, the process has been used extensively for gutta pcrcha, rubber, 
cellulose nitrate and casein (Fisher, 1954). Early progress was mainly based on ‘cut 
and try’ basis and the fundamental process was not fully understood. 

In leccnt years cxpeiimental and numerical studies have been done to under- 
stand the transport phenomena in the extrusion process, resulting in much scientific 
knowledge, primarily for single screw extruders. The first theoretical research on 
plasticating extruder was done by Rowell and Finlayson (1922). Most of the initial 
research (Tadmor and Gogos, 1979, Fenner, 1979, Tadmor and Klien, 1970)weis 
preliminary in nature and relates to the extrusion of polymers. Griffith (1962) 
solved for fully developed flow of an incompressible power-law fluid in a screw 
extruder. The velocity and temperatuie arc essentially the same as those in chan- 
nel of infinite width and length. The effects of curvature and leakage across the 
flights were also ignored. Marshall et al. (1965) and Palit(1974) reported detailed 
experimental results obtained Irom extruders processing materials such as LDPE, 
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polypropylene, polyviuylchloiide, whieh showed screw surface temperatures rising 
rapidly along the feed and compression section to steady values slightly higher than 
the barrel temperature along most of the metering section. Zamodits and Pearson 
(1969) obtained numerical solutions for a fully developed 2-D non-Newtonian flow 
of polymer melts in infinitely wide rectangular screw channels, taking into account 
the effect of transverse flow and superimposed steady temperature profile. Bigg and 
Middleman (1974) and Lidor and Tadmor (1976) have performed theoretical analysis 
to deteiinine the residence time distribution function and strain distribution in screw 
extruders Tadmor and Gogos (1979) and Fenner (1979) have solved the flow of 
a polymer in the feed, compression and metering section of an extruder. Fenner 
(1977) also solved the case of temperature profile developing along the length of 
screw channel. Agur and Vlachopolous (1982) have studied the flow of polymeric 
mateiials, which included a model for the flow of solids in hopper, a model for 
solid conveying zoik' and a mode! for melt conveying zone. Mokhtarin and Erwin 
(1982) developed a mathematical model for mixing in a single screw extruder. Lawal 
and Kalyon (1993) included wall slip in numerical calculation. Gupta et al. (1992) 
have developed a 3-D finite element model for incompressible flow of non-Newtonian 
fluids which can be applied to simulation under isothermal condition. Karwe and 
Jaluria (1990) have numerically analyzed by finite difference techniques the flow of 
polymeric material in metering section of single screw extruder. Sastrohartono et 
al. (1995) dcvclopcxl a numciical model of 3-D transport associated with plastic 
extrusion in a single screw extruder using finite element scheme and compared the 
mimcrical results sa.tisf<ictot ily with experiments done by them. As predicted by 
Das and Ghoshdastidar (1994a,b) the results of Sastrohartono et al. (1995) clearly 
indicate the inadequacies of 2-D modelling. 

The quasi 3-D steady state model developed by Das and Ghoshdastidar (1994b) 
considers cross thermal convection, i.e. convections normal to screw flights and 
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screw roots, and coupled heat tiaiisfer in the melt and the screw. The model is for 
the metering section. The flow is assumed to be hydrodynamically developed but 
theimally undeveloped. 

In food extrusion rheology, the key issue is obtaining the useful data. As pointed 
out by Schwartzberg (1979), the essential difficulty in measuring dough rheology 
is that rheological properties are affected by kinetic phenomenon, such as cross 
linking and protein breakdown. One possible solution, as suggested by Ilohner 
(1978) is to use multiple dies on an extruder. In this case material with identical 
histories would be sub jected to various shear stresses and corresponding shear rates 
could be measured, however its use has not yet been reported. Chen et al.(1978), 
Fong (1978) and Crevone and Harper (1978) obtained flow curves under various 
temperature and flow rate conditions. From flow curves, rheological properties are 
inferred. Harper, Rhodes and Wanninger (1971) suggested a logarithmic mixing 
role form to account for changes in apparent viscosity with changes in moisture and 
the Arrhenius equation form to correlate with changes in temperature. This model 
did not take into account changes in apparent viscosity due to cooking of the food. 
To model time temperature effects on the apparent viscosity of soy protein doughs, 
Remsen and Claik (1978) applied the work of Roller (1979). This model appeared 
useful for initial stages of denaturation only. Morgan et al. (1980) have developed 
a model which <lescril)('s the effects of temperature-time history, temperature, shear 
rate, and moisture content on the apparent viscosity of defatted soy flour dough. 
'Po evaluate this model, temperature as a function of time is required. No single 
model exists desciibing apparent viscosity as a function of composition, variability 
of ingredients, time, temperature and other extrusion parameters. The power-law 
model has been used to describe the isothermal flow through the metering section 
by llarmann and Haiper (1974) and Tsao et al. (1978), and for .cooking extrusion 
by Fricke et al. (1977) Once a simple model that describes the important changes 



in viHcosil.y that. ()((iu wil.li Uinc' (lining (,o(;king of loud (Icjugli oxihUs, it, will be 
tlu'orc'tic'ally possible' to simulate; the; e'litiio cooking extrusion pro(x;ss leading to an 
iiH'u'ase'd uiiele'i sl.anding of the' pimei|)al variabl(;s and their infliie'rie.e on extrusicjii 
rate, cneigy input and finish(;d product characteristics (Harper, 1978). 

1.4 Objectives of the Present Study 

'I'lie objectives of j)i('serit study are twofold. 

First Part 

1 To extend the e^uasi 3-D model of Das and Ghoshdastidar (1994b) to non- 
|)ower law polymer like Vistasil and compare the numerical results with the 
experimental results of Sastrohartono et al. (1995). It may be noted that Das 
and Ghoshdastidar (1991b) did not experimentally validate their numerical 
inode'l llowevei reeently, Das and Ghoshdastidar (1996) successfully validated 
their results for power-law fluids such as LDPE with the experimental work of 
Palit(1974). 

2. 'I'o extend the ejiiasi 3-1) model of Das and Ghoshdastidar (1991b) to a New- 
tonian lluid (which is also a food) such as heavy corn syrup for which the 
viscosity model is available in Sastrohartono et al. (1995) and then compare 
the results lor the same screw configuration and input data (dimensional) with 
LDPE. 

Second Part 

'I'his jiart contains the food extrusion modelling. For this, cooking of defatted soy 
flour dough has been modelled. In the absence of a suitable simple model of viscosity 
for defatted soy flour dough as diseiissed in art. 1.3, it is assumed that through out 
the eooking process the soy Hour dough acts like a polymer lollowing power-law 
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viscosity inoclcl. 'I'lius the quasi 3-1) iiioclcl of Das and Choshdastidar (iOyib) can 
be used. The input data for 25%, 28% and 33% moisture content defatted soy flour 
dough have been taken from Fong (1978) to see how good or how bad is the melting 
a[)ptoximatiou through out the cooking process. While the earlier researchers used 
simple models assuming isothermal flow of power-law food, the present study uses a 
ligorous quasi 3-D model describing non-isothermal flow of food following power-law 
viscosity model. The model also takes into account screw c:onduction 

It may bc' notc'cl that mass tiansfer is not occurring duiing extiusion ol the soy 
flour dough 'Fliis is bc'cause of the fact that the discharge picjssure m food extruders 
ty()ic ally vaiic's bc'twcv'u 30 to (iO atm. At these elevated piessuies, boiling or flashing 
of moisture does not occur within the confines of the barrel because the pressure 
cxcec'ds the vapour [)ic\ssurc of water at the extrusion temjieratc. Once the food 
emerges from the die, the pressuic is released causing the product to expand with 
l.hc' c'xl.c'iisivc' Hashing of moistuie. d'he loss of moisture fioni the jn'odnct lesults 
in adiabatic cooling of the food materials with the product reaching a temperature 
of approximately 80°O in a matter of seconds, where it solidifies and sets often 
retaining its expanded shape (Harper, 1978). 




Figure 1 1: A typical singlr-scrcw extruder 










Chapter 2 


THEORY AND BACKGROUND 

2.1 Screw Extrusion 

In order to predict the flow rates, power consumption, pressure drop, and tempera- 
ture distribution inside a food extruder, the flow pattern or the velocity profile must 
be determined. 

In a screw extruder materials are constantly being conveyed forward through the 
rotation of a close fitting screw inside a barrel. It is the drag of the barrel wall on the 
fluid that will advance the material down the channel in a direction opposite to the 
rotation. This can be visualised by unwinding the screw channel as shown in Fig.2.1 
so that it becomes a straight channel with barrel surface as infinite plane. For ease 
of analysis the coordinate is fixed to the screw and thus barrel moves in a direction 
opposite to the screw rotation and its speed Vj is resolved into two components Hx 
and Vbz, the former along the direction normal to the screw flights and the latter 
along the length of the channel. 

The flow pattern inside the channel is composed mainly of four types of flow. The 
drag flow, is created by the velocity along the helix. It is this component which 
largely drives the fluid forward. The velocity component Vbx creates the transverse 
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flow. It contributes nothing to the pumping rate but will affect the flow pattern and 
is responsible for useful mixing. The die or restriction set up at the exit port of the 
extruder will create the third kind of flow, which sometimes is called back flow or 
pressure flow. For a very small die opening this flow can be in the reverse direction 
and might improve mixing in the system. The last type of flow is the leakage flow, 
which is flow of fluid over the clearance between flights and barrel surface. This flow 
is generally neglected. 

Since our main objective is to predict the pressure and temperature rise which 
occurs primarily in the metering zone of the screw which is also responsible for high 
mixing because of its small depth, in this work the flow and heat transfer in the 
metering section have been analysed in detail. However, a complete model of an 
extruder must include compression as well as feed section. 


2.2 Problem Formulation 

2.2.1 Physical Description of the Model 

The primary difference between a quasi 2-D model and quasi 3-D model is the 
inclusion in the latter of the cross-convection terms i.e., thermal convections normal 
to the screw flights and the base of the screw channel. Similar to the quasi 2- 
D analysis, the flow is assumed to be hydrodynamically developed but thermally 
undeveloped. Both the x-boundaries are assumed to be insulated considering that 
a negligible amount of heat is conducted between the melt and the uncooled screw 
(Fenner, 1977). Also, 2-D conduction is assumed to be present in the x-y plane of the 
screw body with both boundaries in the x-direction acting as insulated considering 
that the screw body is sufficiently thick in that direction. The conjugate heat 
transfer at the melt-screw interface is modelled by assuming the screw body behaving 
like one in which small heat conduction is present in a very thin layer below the 
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scicw suir<u'c while llic rest, ol Uic scicw leinains adiabatic. This assuiription is 
realistic taking into account similar experimental observations by Martin (1970). 

2.2.2 Major Assumptions 

1. h\)i simplilyiiig the analysis, the coordinate system is lixed to the screw and 
thus the barrel moves in a direction opposite to the screw rotation. 

2. The curvature effects have been neglected a.s the screw diameter is large as 
compared to the channel height and therefore the screw is treated as unwound 
(Fig.2.1). 

3. The fluid is consuleied to be incompressible and flow and heat transfer steady 
<uid ((uasi 3-1). 

4. Cood thcimal contact between melt and metal surface exists. 

5 All fluid properties except viscosity are constant. 

6. As the screw channel is shallow and long (i.e. small H/W ratio) and also 
convective effects are small in comparison to viscous effects, creeping flow 
approximations ate made lor the conservation of momentum (Schlichting, 
1979). 

7. The leakage across the sciew flights has been neglected. 

2.2.3 Governing Equations And Boundary Conditions 

Using tlu' above assumptions tin' governing cejuations in the non-dimonsional 
form in terms of the dimensionless variables (listed in the Nomenclature) are 
given below. Tlu' length dimensions have been normalised by dividing them 
by 11, velocities by and pressure by p. 
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Melt 

Continuity: 


du* dv* 
dx* dy* 


X-momentum: 


( 2 . 1 ) 


dx* dx* ^ dx* dy* ^ dy* dx* 


( 2 . 2 ) 


Y-momentum: 


5j/* dy* ^ dy* 5a;* ^ dy* dx* 


(2.3) 


Z-Momentum 




dw* 

dz* 



dw* . 


(2.4) 


Energy: 


n / 

5? + " a? + = 


d^9f d^9f ^ * .,2 r \ 

(2.5) 


Screw 


Energy: 


d^s d^9s 
dx*^ dy*^ 


( 2 . 6 ) 


The viscosity models used in the present study are given in Eq. (2.7), Eq. 
(2.8) and Eq. (2.9). 
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The flow in the screw channel is solved in terms of the cross-channel and 
down-channel flows. The cross-channel flow is defined by u and v velocity 
components in the x-y plane and while down-channel flow is defined by w 
velocity component in the z-direction. The cross- and down-channel flow 
are coupled through y* present in Eq. (2.1)- (2.4) for a non-Newtonian case. 
Here y* is a function of local shear rate 7 *, which is invariant of the rate of 
deformation tensor. In the present case /z* is: 

/z* = /z//z 

where /z and p. in the power law model case are of the form: 

y = \xp[-b{T - To)] (2.7) 

7o 

= ‘==pI-I>{T,-T,)] (2.8) 

in the case of Viscasil are : 


^ 1 + C[Aexp{§)jf~"' 

Aexp{Y) 

1 + C'[Aexp(^)^] 
and in the case of Corn syrup : 


/z = /z<,exp[- 6 (T - To)] 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


p = potxp[-h{T - To)] 


( 2 . 12 ) 
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More details of these viscosity models are given in chapters 3 and 4. 

The boundary conditions are also obtained in the non-dimensional form as: 
at z* = 0, 

u* = u* = = 0 (2.13) 


For 2 * > 0 


at y* = 0, 



at the interface, i.e. 
at y* = 


u* = 0, u* = 0, to* = 0 
‘ ■’ dji- K, dr 


at the barrel surface, i.e. 



u* = v^,o* = 0,to* = 1,^/ = 1 
Vbz 


at X* = 0 and x* = ^ 


(2.14) 


(2.15) 

(2.16) 


(2.17) 


u* =0,u* =0,10* = 0, 


(2.18) 


^ = .^ = 0 

dx* dx* 

The constraint on the flow in the dimensionless form is, 


(2.19) 


rW/H n+Ha/H 


Q 


rW/H fl 

I I w*dx*dy* 

Jo Jh,/h ^ HWVb, 


<lv 


( 2 . 20 ) 
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Here the parameter represents the dimensionless volumetric flow rate also 
called the throughput emerging from the extruder. These boundary conditions 
in the dimensional form are shown as in Fig.2.2. 

Thus for a given extruder, the parameters that govern the solution are: Pe,G,n, 9 „ 
the dimensionless viscosity parameter, [= b(r, — To)] and KslKj. Further- 
more the pressure P*av in Eq. (2.4) is a form of space-averaged pressure over 
x-y plane at any z location. 

2.2.4 Method of Solution 

A finite volume solution methodology has been adopted. The overall method 
of solution of the flow field and temperature has been developed following the 
procedures of Patankar and Spalding (1972), Raithby and Schneider (1979) 
and Vandoormal and Raithby (1984). A fully implicit scheme is used to march 
in z* direction. 

2.2.5 Handling of Melt-Screw interface 

A discretized energy equation is obtained at the melt-screw interface satis- 
fying the energy equation of the melt and screw (Eq. (2.5) and Eq. (2.6) 
respectively) and compatibility conditions (Eq. 2.16) following the method of 
Carnahan et al. (1969). The details can be found in Das (1993). 


2.2.6 Overcdl Solution Algorithm 

The overall solution algorithm is shown by a flow chart given in Fig.2.3. 
The computer program is capable of also predicting correctly the back flow 
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situations due to smaller die openings (low qy cases) that frequently arise in 
case of processing of food doughs. 





START 

V y 



Compute u*,v*,w* & 6p*^/5z* for6^= 0 by solving eq. (1) - (4). 


Guess u*,v*,w*,p*, 0, 5p* /5z* at the next station along z*. 

av 



Solve eq. (2) & (3) for u*,v*using SIMPLE algo, of Patankar (1980). 
Modify pressure u.sing pressure correction term. Correct u*, v*. 

Solve eq.(4) for w*. Correct 5p*/5z*& w* by satisfying eq.(14). 
Solve for0 using first order upwinding in eq.(5). 

Err= max.among diff. of prev. & present value of u*,v*.,p*0. 



Printu*,v*,w*,p,5p* /5z*,6. 

av 













Chapter 3 


VALIDATION OF THE 
NUMERICAL RESULTS 
WITH EXPERIMENTAL 
RESULTS OF 
SASTROHARTONO et al. 
(1995) FOR VISCASIL 


3.1 The Viscosity Model 

In this chapter, the results based on the quasi 3-D model of Das and Ghosh- 
dastidar (1994b) are presented for the fluid Viscasil. Like LDPE , Viscasil is 
also a non-Newtonian fluid. The viscosity model given in Eq.(3.1) for Viscasil 
300-M is taken from Sastrohartono et al. (1995). It is a non power-law model. 
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Aexpi^) 

/• = ( 3 . 1 ) 

1 + C[Aexjp{^)i\ 

where p, is the viscosity in poise, T is the temperature in K, and 7 is the shear 
rate in s The values of constants for the present case are: 

A = 0.7249023 poise. 

B = 2560.804 K. 

C= 7.42082 X 10-® 

power law index n = 0.2671 

Eq.(3.1) is valid over the range of temperatures and shear rates considered 
here. 


3.2 Screw Configuration and Input Data 

The Screw Configuration and input data used in present calculation is given 
in Table 3.1 and the involved dimensionless parameters are given in Table 3.2. 


3.3 Results and Discussion 

The flow and temperature field over the cross-section at any down-channel 
position can be studied. Fig.3.1 and Fig. 3.2 show respectively the develop- 
ment of melt temperature and velocity w* in z direction along the down 
channel direction for three z* locations. It is important to note that 6 and w* 
are the integrated average quantities over x*. Due to the heat input at the 
barrel and the viscous heat dissipation the fluid temperature rises along the 
channel. This results in decrease in the local viscosity ajid hence the velocity 
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Table 3 1. Screw configuration and input data 


Screw Diameter, D 

30.85 mm. 

Pitch 

28.0 mm. 

Helix angle, <f) 

16.1888° 

Maximum channel depth, H 

4.7 mm. 

Maximum channel width, W 

11.52 mm. 

Screw tip width 

1.924 mm. 

Clearance between screw tip and barrel 

0.075 mm. 

Length of metering section, L 

120 mm. 

Density, p 

979 Kg m~^ 

Thermal Conductivity fluid, K j 

0.155758 

Thermal Conductivity screw, Ks (for steel) 

45 

Specific Heat, Cp 

1507.16 3Kg-'^K-\ 

Barrel temperature Tb 

OO 

o 

o 

O 

Screw speed N 

60 RPM 

Inlet temperature for the metering section T, 

53.8° C 

Normalized flow rate 

0.32 


Ta ble 3 2- Dimensionless quantities involv ed 
Peclet number, Pe 4144.176 
Griffith number, G 0.012 


field changes, since or throughput is maintained constant in the extruder 
for steady state condition. 

The velocity vector plots in Fig 3.3 show a strong recirculating flow across 
the screw channel which together with down channel velocity results in a 
spiral motion of the material along the channel. As the material flows the 
temperature rises due to viscous dissipation and heat input from the barrel. 
The screw surface temperature and pressure vs. channel length are plotted 
in Fig. 3.4(a) and Fig. 3.4(b) respectively. It is seen that as expected both 
continuously increase along the screw channel. 
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3.4 Experimental Validation for Various Cases 

The Experimental results presented here were obtained by Sastrohartono et 
al. (1995) using a single-screw extruder with a self wiping screw profile. It is 
to be noted that the screw configuration is same as in Table 3.1. 

In the first case barrel temperature Tj = 25.0°C, inlet temperature T, = 18.9°C, 
normalized flow rate = 0.34 and screw rotation speed N = 10 rpm.(See Table 
3.3). The outlet temperature is the bulk temperature at the exit. 

In the second case comparison was done with barrel temperature Tb = 80°C 
(Table 3.4). The results shown in Table 3.3 and Table 3.4 indicate clearly that 
the quasi 3-D model is highly reliable. 


Table 3 3 Comparison between computed model and experimental results for Viscasil, N = 10 
rpm, Tt = 25 0°C, T. = 18 9”C 


Test 

Outlet Pressure 
[bars] 

Outlet 

Temperature ["Cj 

Experimental 

Results 

4.750 

24.20 

Computed 

Results 

4.854 

25.25 
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Table 3 4 Single-screw extruder characteristic experimental and computed results Tj, = SO^C, T, 
= 53 8°C; Viscasil-300M 


Mass flow 
Rate 

IKgft-'l 

Screw speed 
[rpm] 

Outlet Pressure 
(Experimental) 

[x 10®iVm-2] 

Outlet Pressure 
(Computed) 

[x 10®iVm-2] 

1.0 

10 

8.60 

8.94 

3.6 

36 

14.31 

14.62 

6.0 

60 

18.74 

18.96 

9.0 

90 

i 

22.44 

22.59 




z*=45.79 

— T , , 

S 

1 1 

^ 

N 

* 

II 

VO 

in 

00 


t'3 1 Toniperalurc profiles along the screw channel of a single-screw extruder Fluid Viscasil 

80“(1. 7 ; = 53 8"C, = 0 32, N = ODrpin 






F.jurc 3 2 .>• prof.l.= along the see channel of a s,ngle-scre« extrude, Flu.d \', seas, I, % 
80" C\ 7. = 53 8"C, -7„ = 0 32, N = 60rpm 
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Figure 3.4 Variation of (a)Screw surface temperature and (b)Pressure along the screw channel of 
a single-screw extruder. Fluid Viscasil, Th = 80^C, T, = 53 8®C, = 0 32, N = 60rpm 





Chapter 4 


COMPARISON OF 
NUMERICAL RESULTS FOR 
CORN SYRUP AND LDPE 

4.1 The Viscosity Models for Corn Syrup and 
LDPE 

In the present chapter, a Newtonian fluid such as corn syrup is considered. 
The numerical results are compared with those for a non-Newtonian fluid 
LDPE (Low Density Polyethylene) for the identical screw configuration and 
input data in the dimensional form. The viscosity model for corn syrup is 
given 35 (Sastrohartono et al., 1995): , 

fi = 1052.8exp[— 0.095(r — 20)] (4.1) 

where fx is the viscosity in poise, T is the temperature in °C. 
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'I'lu- (U'n.sity, Ukmmi.U c oaducUvity cuid Uic apecilic hc<iL oi com syiup arc : 
p = 1381 kg m~'* 

Kf = 0.317 W7n-‘A'-‘. 

(•V = 2015 

While 1,1 k; Viscosily model (power law) for LDPE is : 

A = Ao(~) cxp[-b{T-Ta)] (4.2) 

7o 

wIk'ic // IS t.li(' visco.sil.y in Pa-s, '1' is 1,1k; l,('nipcral,i]te in "(J, and 7 i.s the shear 
laU' in .s ' . 1 h(‘ values ol constants lor the present case are; 

p„ = 2000 Pa-s. 

7 ; = 200 °C. 

b ^ 0.01 

7„ = 1 .s'h 

pow<‘r law iiuh'X n = 0.48 

and (h'lisity f) -- 979 kg ni"'* 

'I'Ik' tlK'rnial conductivity and the specific heat are : 

Ay 0.155758 Wm-i/W*. 
r.; = 1507.16 Jkg-^K-\ 

4.2 Screw Configuration and Input Data 

'I'Ik’ .screw (ilements dimensions are similar to one shown in Table 3.1. The 
dinu'nsionlesH parametens (Pc and G) for each case arc given in Table 4.1 and 
Tal>l(“ 4.2. Not<' the values are different. 
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I'lie case presented here is for: 
barrel temperature Tt = 80° C 
s<'rew speed = 60 RPM. 

inlet temperature for the metering section T, — 53.8° 


Tabl e 4 1- Dimensionless quantities involved for Corn S yrup 

Peclet number, Pe 3281.265 

Griffith number, G 0.9072 

Dimensionless throughput, 0.32 


Ta ble 4 2 Dimensionless quantities involved for LDP E 
Peclet number, Pe 2734.176 

Griffith number, G 2.0134 

Dimensionless throughput, 0.32 


4.3 Results and Discussion 

i he comparisons of w* velocity and melt temperature at the various down 
channel positions are plotted in Fig.4.1 and Fig.4.2 respectively. From the 
nudt temperature plots it is observed that for heavy Corn syrup the melt 
temperature rise is lower than that for LDPE, which is due to the higher 
viscous dissipation in LDPE as it has higher viscosity. The w* velocity plots 
(Fig.4.2) at the entry and exit of metering section show that at low y* velocity 
for corn syrup is lower than that for LDPE while at higher y* it is just the 
opposite. This can be explained by the fact that for LDPE viscosity is a 
function of temperature and shear rate, while for corn syrup it is a function 
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()1 U'lnpt'i al,m<' t)nly Near the wall the shear rate is higher than farther away 
lioiu the wail. '1 hus lor LDPE, viscosity which is a function of negative power 
ol .sIkmi iat(', is h'sscu' iK'ai the wall than farther from it and the velocity 
((>* ol LDPE will 1)<> highc'r there, while corn syrup has no such effect as it 
IS ui<ie|)<Mi<hMit ol slu'a.i latc'. Due to heat in|)ut at the bariel and viscous 
dissi[)atu>n the fluid temperature rises along the length of the channel. This 
K'sults in d('cr<'a..s(' in local viscosity and hence the velocity field changes. But 
(iie < liangc' ohservi'd is not much since the temperature diflerence is small. The 
s< K'w suiLhc' t('mperatur(' (hdg.d.3(a)) is seen to rise continuously for both 
< a.ses. 11 is int.t'U'sting tt) luite litthi dilference between the two. One reason 
may Ix' tin' almost identical thermal conductivity ratio (K^IKf) for LDPE 
and coin syiup, as lor both cases the same screw material is taken and A/ for 
both aie almost the same. 'Phe slightly higher screw surface temperature for 
coi n syt np in a, bout. Hist (>()% length of the channel may be due to the fact that 
b«'( aus<' of higher cross thermal convection for corn syrup (large Peclet number) 
tin* barrrd h(>at is convected better to the screw surface and hence the screw 


sml'a<-e temperature is greater. However, towards the exit the screw surface 
lempeial.uK' for th<' com syrnp is lower because of greater viscous dissipation 
(>lle< I. for LDPE a.s it has higher aiiparent viscosity (Griffith number for LDPE 
IS guMler). Also, it is noU'd that the screw surface temiierature for both the 
< ases appioaelu's the barrel temperature, due to the combined effects of cross 
thermal eonv<«ction and viscous Imat dissipation. The pressure is seen to rise 
continuously for both cases (Eig.4.3(b)). But the rise is higher for LDPE than 
fot corn synii), again due to higher vi.sc.osity of LDPE. 
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Figure 4 3: Comparison of (a) Screw surface temperature and (b) Pressure along the screw channel 
between LDPE and corn syrup 





Chapter 5 


NUMERICAL RESULTS FOR 
DEFATTED SOY FLOUR 
PROCESSING AND ITS 
VALIDATION WITH 
EXPERIMENTAL RESULTS 
OF FONG(1978) 


5.1 Viscosity Model used for Defatted Soy Flour 
Dough 


In this chapter defatted soy flour is considered and a comparison of the results 
is done with the experimental results obtained by Fong (1978). 
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The viscosity model used for defatted soy flour is power law model, given as 
follows: 


= ^exp[-h{T-To)] (5.1) 

Ho 

where (i is the viscosity in Pa-s, T is the temperature in ^C, and 7 is the shear 
rate in 

The power-law index, n and constant, calculated from the experimental 
results obtained by Fong (1978) are listed in Table 5.1 for three moisture 
contents. The method used for computation is given in Appendix A. 


Table 5 1 : The value of /Xp and n for three moisture contents 


Dough flour 
moisture 
[% by weight] 

33% 

28% 

25% 

1^0 

[Pa-s] 

32763.0 

123010.0 

134215.0 

n 

0.4569 

0.1222 

0.0177 


For each case, the value of b and % are O.OVC ^ and Is ^ respectively. 


5.2 Screw Configuration and Input Data 


The results presented here for quasi 3-D model are compared with the exper- 
imental results obtained by Fong (1978). The screw configuration and input 
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data are given in Table 5.2. While the dimensionless parameters in three cases 
are given in Tables 5.3, 5.4 and 5.5. It may be noted that in absence of reliable 
information regarding ‘b’ for soy flour , the corresponding value for LDPE has 
been used. Similarly, the thermal conductivity, Kj for soy flour dough was not 
available and hence Kj for wheat flour (Singh, 1978) was used. It is expected 
that aforementioned uncertainties in input data will not seriously effect the 
results as for the soy flour processing ‘melting’ model is used and conductivity 
of wheat flour should be close to that of soy flour. 


Table 5 2 Screw configuration and input data 


Screw Diameter, D 

25.4 mm. 

Helix angle, (j) 

17.5° 

Maximum channel depth, H 

1.46 mm. 

Maximum channel width, W 

19.3 mm. 

Length of metering section, Lm 

167.64 mm. 

Density, p 

1220 Kg m-2 

Thermal Conductivity (dough), Kj 

0.45 Wm-^K-L 

Thermal Conductivity (screw), K, 

45 Wm-^ir-T 

Specific Heat, Cp 

3320 iKg-^K-'^. 

Barrel temperature Th 

73.89° C 

Screw speed N 

26 RPM. 

Inlet temperature for the metering section T) 

62.7° C 


Table 5. 3’ Dimensio nless parameters involved for 33% moistu re content Soy flour dough 
Peclet number, Pe 396.260 

Griflflth number, G 1.588 

Dimensionless throughput, 0.0942 
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'Table 5 4 Dimensio nless parameters involved for 28% moisture conte nt Soy flour dough 

Peclet number, Pe 396.260 

(iiillith luunbcr, G 2.0393 

Dimensionless throughput, 0.1847 


'Table 5 5 Dimensio nless parameters involved for 25% moisture conte nt Soy flour dough 
IV'clct iiumbci, I’c 396.260 

Grillith muiiber, G 1.5915 

Dimensionless throughput, qy 0.1955 


5.3 Results and Discussion 


'I'lie velocity vector plots for 25% moisture content dough (Fig. 5.1) show weak 
lecirculation in the cross-sectional plane which is due to the high viscosity of 
low moisture dough. A high recirculation is seen for 28% and 33% moisture 
content dough ( 1' ig. 5.5 and Fig .5.9 respectively) since viscosity is compar- 
atively lower in this case due to high moisture presence which facilitates the 
motion. I'he down channel velocity w* variation with y* at various z locations 
shows that back How occurs. This was also reported by Fong (1978). The 
baddlows lor 25%(Fig. 5 2) case and 28% (Fig. 5.6) case are small while for 
the case of 33% moisture backflow is larger. This is due to the difference in 
normalised flow rate of the three cases (see Table 5.3, 5.4 and 5.5). Since a 
smaller </„ imiilies that the die ojieiiing is small, thus not all dough can pass 
through it and some of it flows in the backward direction. 

tl.e temiieiatuK' prolih- at vaiious z locations, it is noticed that temper- 
ature is always higher than the barrel temperature. This was also reported 
during the experiments by Fong (1978). Also temperature rise is higher in the 
case of 25% dough (Fig. 5.3) than other two cases (Fig. 5.7 and Fig. 5.11). 
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This IS due lo liigliof moisture contents and lienee lesser viscous dissipation in 
the latter cases. 

'I'lu' picssuie is .sec'ii to rise linearly for all three cases viz. 25% (Fig 5.4(b)), 
28% (h'lg. 5 8(b)) and 33% (Fig. 5.12(b)). Though pressure rise is more in the 
t ase ol 25% moisture than 28% case, which in turn is more than the 33% case 
due to higher viscosity for the former cases. The screw surface temperature 
list's vt'ry rapidly in the early part of the metering section and then steadies 
to a temperature much greater than the barrel temperature (Fig. 5.4(a), Fig. 

5. 8(a) and Fig. 5.12(a)). 'I'his is due to very high viscous heat dissipation in 
food as compared to polymer 

5.4 Experimental Validation 

'Table 5.6 shows computed and experimental outlet bulk temperatures for all 
three cases. It is seen that for low moisture percentages the difference in 
temperature is much higher than that at high moisture percentage. A similar 
tiend IS observed m the case of outlet pressure (comparisons shown in Table 

5.7). 

'I’his can be explained due to the fact that at lower moisture percentage 
viscosity is high and hence more heating takes place which initiates cooking 
reaction. 'This cooking reaction changes the viscosity which can not be simply 
represented by the power law model. But for high moisture percentage case a 
better agi cement bi'tween two results is seen i.e. cooking here has less effect 
on the viscosity of the fluid . So this model is inadequate for high temperature 
and low moisture case but can be effectively used for low temperature and 
high moisture ca.se. 
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Table 5.6 Comparison of computed and experimental outlet temperature 



Table 5 7 Comparison of computed and experimental outlet pressure 


Moisture 

Outlet Pressure 

Outlet Pressure 

percentage 

(Experimental) 

(Computed) 


[x 10® Nm-2] 

[x 10®Nm"^] 

25 

5.040 

2.76 

28 

4.912 

2.67 

33 

1.480 

1.34 
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Figure 5 4 (a)Screw surface temperature and (b)Pressure variation along the length of screw 
extruder for 25% moisture content dough 








Figure 5 6 le' velocity profiles for 28% Jnoisture content dough at three down channel locations 
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Figure 5 7 lernperaturc profiles for 28% moisture content dough at three down channel locations 
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Figure 5 8 (a)Screw surface temperature and (b)Pressure variation along the length of screw 
extruder for 28% moisture content dough 









Figuro 5 10 IV* velocity profiles for 33% moisture content dough at three down channel locations 
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Figure 5.12: (a)Screw surface temperature and (b)Presaure variation along the length of screw 
extruder for 33% moisture content dough 





Chapter 6 

CONCLUSIONS AND SCOPE 
FOR FUTURE WORK 

6.1 Conclusions 

A numerical study of 3-D flow and heat transfer during processing of polymers 
and foods in the metering section of a single screw extruder is carried out. 
A quasi 3-D finite volume modelHng is used and flow field in cross-sectional 
plane is calculated using SIMPLE algorithm. A fully implicit scheme is used 
to march along the length of the channel. The computer model is capable of 
predicting the flow even for very small die opening (very low cases) which 
gives rise to back flow situations, as in processing of food doughs. 

The polymers considered are LDPE and Viscasil, while the foods are heavy 
corn syrup and defatted soy flour. One of the achievements of this work is to 
successfully extend the quasi 3-D model of Das and Ghoshdastidar (1994b) for 
power-law fluids to non-power law fluids like Viscasil. This has also been 
validated by comparison with experimental results of Sastrohartono et al. 
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(c) lb liud that powor-law viscosity model (“melting” model) is reliable 
I'or high moisture content defatted soy flour dough but not so for drier 
doughs. It may also be noted that for the first time a quasi 3-D model 
using powci-law viscosity model has been used to simulate food dough 
piucessmg. 


6.2 Scope for Future Work 

l‘\)r the purpose of initial application, the results have been shown for 
only one dough material (i.e. defatted soy flour) with only three levels of 
moisture content. Once a good viscosity model for all levels of moisture 
contents is developed, it is desirable to extend the applications to other 
food materials such as corn, wheat, rice etc. 



Appendix A 

PROCEDURE FOR 
COMPUTATION OF 
CONSTANTS OF POWER 
LAW MODEL 


A.l Procedure for Calculation of fXo and n 


Here a calculation procedure proposed by Rogers (1970) is presented to 
estimate Ho and n, the constants of power law equation for describing the 
viscosity behaviour of defatted soy flour. The use of procedure requires 
pressure drop and flow data using a number of dies having varying L/R 

ratios. 

The viscosity model used for deflated soy flour is power law model, given 
as follows: 
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M = //„(-)” WH(r-r„)l (A.i) 

lo 

where // is the viscosity in Pa-s, T is the temperature in °C, and 7 is the 
shear rate in 

The experimental data used here for calculation is taken from Fong 
(1978). The dimension of dies used by him for extrusion were: 



The data taken for dough is tabulated in Table A.2 and Table A. 3 for die 
1 and die 2 respectively. 
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For the solution from the data we calculate 7a given as: 


• _i^ 

ttB? 


(A.2) 


The calculated values are given in Table A.4.Then log7o vs. logp is 
plotted and from the graph AP is found for 7a = 1^"^ and for each 
die. These values are given in table A.5 along with L/R values of both 

the dies. 
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Table A 5 Value of AP for -ya Is ^ and 2s~^ for Die 1 and Die 2 


7a 

log AP 
L/R = 9.75 

log AP 

L/R = 9.78 

1 

14.2556 

13.9728 

2 

14.6031 

14.6104 


Now AP vs. L/R is plotted and then extrapolated to obtain I/*/R, which 
is the value corresponding to AP = 0 (given in Table A.6). 

Table A 6' Z.*/R value for 7 a l.s~^ and 


7a 

[^-^] 

L*IR 

1 

8.5192 

2 

9.0925 


These L*/R values are plotted against % to get their values at data 
points. These values are given in Table A. 7. 

The values so obtained are used to calculate t^, which is given as: 

^ (A.3) 

“ ~ [2(s + m 

These values are given in Table A.8. 

Slope of logr^ vs log 7a gives the flow behaviour index. This is found 
out to be n = 0.4569 in this case. For calculation of fio we get — 
(3n+l)(4ii)7a and which is used to get fi = (given in Table A.9. 
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Tabic A 7- L* /R value for data points for Die 1 and Die 2 


Reading 

L*/K 

Diel 

L*IR 
Die 2 

1 

8.4497 

8.4176 

2 

8.9858 

8.7026 

3 ^ 

9.3490 

9.0544 


Reading 

Ty, [xlO^ Pa] 
Diel 

Tyj [xlO^ Pa] 
Die 2 1 

1 

4.0669 

2.8000 

2 

5.5771 

4.0660 

3 

6.6659 

5.5186 


The intercept of log/i vs. log 7 gives the value of fJ-o- 
calculation /io 3.2763 x 10 Pa. 


In the present 
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1 and Die 2 
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